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A weakly random Universe? 
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ABSTRACT 



Cosmic microwave background (CMB) radiation is characterized by well-established scales, the 2.7 K temperature of 
the Planckian spectrum and the 10 amplitude of the temperature anisotropy. These features were instrumental in 
indicating the hot and equilibrium phases of the early history of the Universe and its large-scale isotropy, respectively. 
We now reveal one more intrinsic scale in CMB properties. We introduce a method developed originally by Kolmogorov, 
which quantifies a degree of randomness (chaos) in a set of numbers, such as measurements of the CMB temperature 
in a given region. Considering CMB as a composition of random and regular signals, we solve the inverse problem of 
recovering of their mutual fractions from the temperature sky maps. Deriving the empirical Kolmogorov's function in 
the Wilkinson Microwave Anisotropy Probe's maps, we obtain the fraction of the random signal to be about 20 per 
cent; i.e., the cosmological sky is a weakly random one. The paper is dedicated to the memory of Vladimir Arnold 
(1937-2010). 
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The high-accuracy Planckian spectrum, the quadrupole 
anisotropy, and the acoustic pea k s of the power spectrum 
(|Mather 19941 ISmoot et al"l992l Ide Bernardis et al 2000|) 
of the cosmic microwave background (CMB) appear to 
be carriers of valuable information on the early Universe. 
Among other peculiarities of the CMB signal is that it is 
described well by a Gaussian distribution of initial fluctua- 
tions. This is remarkable at least for two reasons. Gaussian 
fluctuations follow from the simplest single scalar field ver- 
sion of inflation. The Gaussian is also a limiting distribu- 
tion, in accordance with the central limit theorem, as the 
sum of many independent random variables of finite vari- 
ance. CMB fitting the Gaussian certainly does not reveal 
its randomness, as it might arise from both random and 
regular distributions. 

Is it possible to reveal the fractions of random and reg- 
ular signals in the CMB? 

One can deal with this problem using the Kolmogorov 
distribution and the stochasticity parameter as a rigorous 
measure of randomness (jKolmogorov 1933} lArnold 20 08a). 
Arnold applied this technique to arithmetical and geo- 
metrical progressions and to number sequences (Arnold 
2008abc,2009ab), and it has also been applied to CMB tem- 
perature datasets ( |Gurzadyan et al 2009[ ). Now, based on 
the 7-year CMB data obtained by the Wilkinson Microwave 
Anisotropy Probe (WMAP) (jJarosik et al 2010|) we have 
derived a general result, i.e., that the random signal yields 
about 0.2 fraction in the CMB total signal. 

We start by briefly outlining the approach. We consider 
a sequence {Xi, X2, ■ ■ ■ ,X n } of a random variable X or- 



dered in an increasing manner X\ < X2 < • • • < X n . 
Cumulative distribution function of X is defined as the 
probability of the event X < x 

F(x) = P{X < x}, 

and an empirical distribution function F n (x) as 

F n (x) = , forx < X x ; 

k/n , for X k < x < X k+ i, k = 1, 2, . . . , n - 1 ; 
1 , for X n < x . 

Then the Kolmogorov's stochasticity parameter A„, which 
is also a random quantity, is defined as the normalized de- 
viation of those two distribution functions 

A„ = V^ sup\F n (x)-F{x)\ . (1) 

X 

Here, the multiplier y/n appears to account for the devia- 
tion expected for n independent observations of a genuinely 
random variable. 

The remarkable result is Kolmogorov's (1933) proof 
("the astonishing Kolmogorov's theorem" as Arnold 
(2008b) called it), concerning the asymptotic behavior of 
the distribution of A for any continuous F 

Km P{X n < A} = $(A) , 

n— J-oo 

where the limiting distribution for any real A > is 

+00 

$(A) = (-l) fe e _2fe2A2 , *(0) =0, (2) 
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the convergence is uniform, and <E>, the Kolmogorov's dis- 
tribution, is independent of the distribution F of the initial 
random variable X. The interval of probable values of A 
yields 0.3 < A„ < 2.4 (Kolmo gorov 1933] ). This universal- 
ity of Kolmogorov's distribution marks it as a measure of a 
stochasticity degree of datasets (Arnold 2008abc,2009ab). 

This objective measure enables one to consider 
the composition of signals of various randomness; e.g., 
the behavior of A„ and $ was studied at numeri- 
cal simulations of sequence s ( Ghahra manyan et al 2 009 , 
|Mirzoyan fc Poghosian 2009[ ). 

Within this remarkable system, an approach to ap- 
ply is the CMB dataset, just because its cumulative dis- 
tribution function happens to be known, and it is close 
to Gaussian! Then for Gaussian distribution function F 
and via a properly developed strategy (for algorithmic 
and numerical details see Gurzadyan and Kocharyan 2008; 
Gurzadyan et al 2009) one can obtain the distribution 
of 3> for CMB, the Kolmogorov map for WMAP's tem- 
perature data, as shown in Fig.l. Not only is the non- 
cosmological signal, the Galactic disk, clearly outlined, 
which was not at all trivial a priori, but also a non-Gaussian 
feature such as the Cold Spot also appears to be out- 
lined by a specific behavior by the function $ that is com- 
patible with its void nature (Gurzadyan and Kocharyan 
2008,2009; Gurzadyan et al 2005,2009). The Gaussian na- 
ture of CMB is characterized by certain correlations stud- 
ied by traditional methods, i.e. the harmonic analysis, cor- 
relations functions, and power spectrum, which were par- 
ticularly informative due to the structures of the acous- 
tic peaks. Application of Kolmogorov's approach will com- 
plement those studies by quantifying the cumulative de- 
gree of correlations in the CMB signal. Kolmogorov's pa- 
rameter is also efficient in detecting point sources, radio 
and Fermi/LAT gamma-ray, and other anomalies in CMB 
maps ( Gurzadyan et al 2010| [Gurzady an et al 2008[ ); i.e., 
this method enables the separation of the cosmological sig- 
nal from non-cosmological ones. 

In the present paper, for the first time, we undertake a 
novel analysis of the cosmological CMB signal revealed by 
the $ distribution. First, we obtain the Kolmogorov dis- 
tribution $(A) for different regions of sky and WMAP's 
7-year W-ba nd data (| Jarosik et al 2010|) in HEALPix rep- 
resentation (IGorski et al 2005]) . The procedure was as fol- 
lows. Circles of radii 1°.5 were considered to have enough 
pixels to estimate the mean $(A), given that the A it- 
self is a random parameter. The variance of the Gaussian 
then was estimated for each circle region containing around 
on average 540 pixels. Then the resulting frequency dis- 
tribution vs $ was obtained, and it appears to be decay- 
ing, corresponding to 10000 2D-balls randomly distributed 
over the sky with the consequent elimination of belts of 
coordinates |6| < 20°, \b\ < 30°, \b\ < 40°, in order 
to evaluate the role of the contamination of the Galactic 
disk. Figure 2 shows the plot for \b\ < 30°. The 2D- 
ball s obviously also cover the cold spot reg ions studied 
in (jJarosik et al 20101 |Gurzadyan et al 2009|). The p oint 
sources revealed in CMB maps ( |Gurzadyan et al 2010| ) are 
also covered; however, they cannot contribute to the distri- 
bution of $, not only because of their small quantity but 
also mainly because of their point-like feature, i.e. small 
number pixelization. 

The first conclusion drawn from Figure 2 is expected. 
The CMB signal is certainly not purely random Gaussian; 




Fig. 1. The Kolmogorov map of cosmic microwave back- 
ground radiation, i.e. the distribution of Kolmogorov's 
function <£> over the sky for the Wilkinson Microwave 
Anisotropy Probe's 7-year W-band temperature dataset. 
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Fig. 2. The empirical frequency distribution of 
Kolmogorov's function <i> in the CMB map with the 
extracted Galactic disk region |6| < 30° and the generated 
signal z n for the fraction of the random sequence a = 0.19. 



otherwise, the stochasticity parameter, hence, the function 
$, would have a uniform distribution, i.e. constant, as fol- 
lows from Kolmogorov theorem. 

Once Kolmogorov theorem enables estimation of the 
randomness $ for a given sequence, with the empirical dis- 
tribution of $ in Figure 2, one can inquire into the inverse 
problem of reconstructing the initial signal to result in that 
distribution. 

To solve this inverse problem, we consider a regular se- 
quence represented in the form (cf. Arnold 2008c) 



(a + hn) mod(6 — a), 



(3) 



i.e. as an arithmetical progression repeated within an inter- 
val [a, 6], and where h is an arbitrary constant. 

Since the sequence (3) corresponds to a uni- 
form distribution, one can use von Neumann's method 
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Fig. 3. The confidence level in % 2 vs the parameter a, i.e. 
the fraction of the random component, for the WMAP7 
CMB maps with extracted Galactic disk regions: \b\ < 
20°, 30°, 40°. 

(jvon N cuman n~T951|) to extract a Gaussian nonrandom se- 
quence of standard deviation a and mean value m from 
it. The latter we consider within the adopted interval 
a = m — lOer and b = m + lOc. 

Von Neumann's method enables generating of sequences 
with any differential distribution f(x) and extracting of se- 
quences with the desired distribution from a uniform one. 
The method is particularly efficient when no explicit form 
of F(x), i.e. of the integral of f(x), is available. The proce- 
dure is as follows. First, we consider a uniform sequence a n 
and a distribution f(x) having minimal and maximal values 
y m in and y m ax- Then, to extract a sub-sequence from a n 
one has to go through all elements of a n and, for each a,, 
generate a random uniform number z within the interval 
(Vmin, Umax); if 2 < /(a»), one assigns Oj to a new sequence 
and skips it otherwise. The extracted subsequence will have 
a distribution f(x). We used this procedure to generate ran- 
dom and regular sequences with Gaussian distribution, i.e. 
when f(x) is a Gaussian. If the original sequence a n is a 
random one, then one gets a random Gaussian sequence, 
and a non-random one otherwise. 

By generating many regular sequences with a Gaussian 
distribution of mean and standard deviations correspond- 
ing to the WMAP's CMB signal and calculating the 
Kolmogorov parameter for each of them, we were able to 
construct the empirical distribution of $ at the best choice 
of the parameter h = 7.342.... 

As the next step, the following sequence was considered, 

z n = ax n + (1 - a)y n , (4) 

where random Gaussian sequence with the same 

m and a as the regular von Neumann-Gaussian sequence 
y„, i.e. the one derived by means of the von Neumann 
method. Thus, we have a sum of random and regular se- 
quences with their varying proportion defined by the pa- 
rameter a € [0, 1]: for a — we have a regular sequence, 
while for a = 1 it is random sequence. 



It can be shown that z n also has a Gaussian distribu- 
tion with mean m and standard deviation yl — 2a + 2a 2 a. 
While varying a from to 1 we see, as expected, the distri- 
bution of $ gradually changing from the form in Figure 2 
to a uniform one. At some value of a, we get the closest fit 
evaluated in confidence levels for distribution of $ of z n to 
WMAP's sequences. Figure 3 shows the behavior of \ 2 y s 
the a for the CMB maps with extracted regions |6| < 20°, 
|6| < 30°, \b\ < 40°, and as seen, the results of the last two 
maps are identical. We see that the minimal values of \ 2 , so 
the highest confidence levels are reached for 0.15 and 0.19 
values of a for the 20° and 30° disks, respectively. 

The resulting a distributions are shown in Figure 2 
along with WMAP's frequency plot for the CMB map with 
extracted Galactic disk of |6| < 30°. If the small difference, 
a few per cent, between the values of a of the CMB maps 
with extracted \b\ < 20° and |6| < 30° regions stems from 
the Galactic disk contamination, then one can state that 
the latter contributes to a slight regularization of the CMB 
signal. 

The methodology of Kolmogorov's approach is based on 
studying statistical properties, namely, the degree of ran- 
domness of the signal via the stochasticity parameter and 
the distribution As it appears, those characteristics are 
informative enough to distinguish one signal from another, 
for example, the cosmological CMB and the Galactic syn- 
chrotron one. Although the photons of such different origins 
can have the same temperature, they differ in their degree of 
randomness, hence in their correlation properties; e.g., the 
synchrotron is characterized by a power-law spectrum, and 
the CMB by a Planckian one. This difference is reflected in 
the statistical analysis. However, here we went further to 
quantify the contribution of the correlations and of chaotic 
components. The correlations in the CMB signal are car- 
riers of information on various processes, from primordial, 
through the last scattering and reionization epochs up to 
the integrated Sachs- Wolfe and Sunyaev-Zeldovich effects. 
Revealing the degree of correlations in CMB signal, partic- 
ularly, reflects the cumulative contribution and properties 
of those processes. 

Thus, considering CMB as a composition of signals, 
we showed that the random component is a minor per- 
turbation of mostly regular signal. The model-independent 
character of this result comes from the universality of 
the Kolmogorov theorem that we used in analysing real 
datasets, the WMAP's 7-year CMB maps. We solved the 
inverse problem by numerical methods and recovered the 
mutual fractions of random and regular signals which would 
lead to the empirical Kolmogorov's distribution. 

That the CMB sky is a weakly random one and the de- 
rived 0.2 fraction of the random signal in CMB by their 
intrinsic statistical content may be compared with other 
scales of CMB, the 2.7K temperature of the Planckian spec- 
trum, and the 10 -5 amplitude of the quadrupole anisotropy, 
indicating the hot equilibrium history of the Universe and 
the isotropy of the Universe, respectively. If the Planckian 
spectrum and the anisotropy amplitude do reflect and do 
scale certain correlations in CMB signal, the obtained frac- 
tion scales the correlations in the cosmological Gaussian 
perturbations. 

What insights will the new scaling open to the early 
Universe? 
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We dedicate this paper to the memory of Vladimir 
Arnold (1937-2010), recalling the memorable discussions 
(VG,AAK) with him. 
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